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P RE FACE 


Hie  research  reported  in  this  Memorandum  is  a  part  of  a  continuing 
investigation  of  the  effects  of  ground  shock  on  underground  structures. 
This  study  was  prompted  in  general  by  a  questioning  of  the  accuracy  of 
free-field  stress  measurements,  and  in  particular  by  an  interest  in 
the  effects  of  gaga  design  on  these  measurements.  Most  of  the  current 
free-field  stress  data  available  are  derived  from  pa  icie-velocity  mea¬ 
surements;  very  few  are  obtained  through  actual  stress  gages,  a  long- 
recognized  inadequacy. 

It.  is  therefore  the  purpose  of  the  Memorandum  to  investigate  the 
feasibility  of  an  or.trM directional  stress  gage  and  to  provide  theoreti¬ 
cal  background  for  the  design  and  use  of  such  a  device. 

Doctor  F.  C,  Moon  is  presently  affiliated  with  Princeton  University 
and  is  also  a  consultant  to  the  RAND  Corporation. 
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SUMMARY 

This  Memorandum  presents  an  analysis  of  the  transient  response  01 
the  pressure  in  an  embedded  elast^  inclusion  due  to  an  incident  com- 
pressicnal  wave.  It  represents  a  continuing  effort  at  RAND  to  evalu¬ 
ate  transducer  diffraction  effects  on  ground  stress-wave  measurements. 
Included  is  an  outline  of  a  design  for  an  omnidirectional  pressure 
transducer  with,  a  response  capable  of  being  interpreted  in  terms  of 
the  theoretical  solution  presented  in  this  Memorandum. 

Both  the  frequency  response  and  transient  behavior  are  treated. 

It  is  found  that  the  pressure  or  mean  stress  at  the  center  of  the  in¬ 
clusion  will  be  insensitive  to  the  curvature  of  the  incident  wave. 

The  primary  source  of  distortion  between  the  inclusion  pressure  and 
the  free-field  pressure  in  the  incident  wave  is  internal  reflections 
in  the  inclusion.  An  early  time  analysis  reveals  that  these  can  be 
minimized  by  matching  the  acoustic  impedance  (product  of  density  and 
coaipressional  wave  speed)  of  the  inclusion  with  that  of  the  matrix. 

An  estimate  of  the  time  for  these  reflections  to  decay  due  to  radia¬ 
tion  "damping11  is  found  to  depend  on  the  impedance  ratio. 

TWo  methods  are  presented  for  obtaining  the  total  pressure  re¬ 
sponse  due  to  a  nonperiodic  incident  wave.  The  fi.  t  uses  the  calcu¬ 
lus  of  residues  and  the  high-frequency  response  to  sum  the  resulting 
infinite  series.  The  second  method  proceeds  by  solving  the  inverse 
problem:  to  find  the  incident  wave  pressure  in  terms  of  the  inclusion 
pressure.  This  latter  method  is  presented  in  the  form  of  an  intenral 
equation  for  the  Inclusion  pressure  which  enables  an  exact  explicit, 
solution  to  be  found  between  successive  reflections. 
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I.  INTRODUCTION' 


/-.mong  the  many  effects  associated  with  a  nuclear-explosion-induced 
ground  ‘hock  that  must  do  considered  in  the  design  of  a  survivable  sys¬ 
tem  is  the  magnitude  of  the  stress  in  the  ground  shock.  In  fact,  most 
design  specifications  of  a  hardened  system  usually  call  for  structures 
to  withstand  a  certain  stress  level.  Thus  the  abilities  to  predict  and 
to  measure  free-field  stress  are  of  prime  importance  in  any  nuclear  ef¬ 
fects  tests. 

Most  of  the  existing  stress  gages  are  devised  primarily  for  static 

(1) 

measurements  and  may  not  be  applied  to  dynamic  stress  measurement . 

There  nave  been  many  recent  studies  attempting  to  produce  better  stress 

gages  capable  of  measuring  high-intensity  ground  stress ,  but  usually 

because  of  their  comp  lex  geometric  shape,  a  precise  analvsts  of  the 

wave- transducer  interaction  is  impossible :  hence,  the  interpolation 

(i  ) 

of  data  Is  in  doubt.  A  strain-sensitive  transducer  of  simple  shape 

and  construction  is  needed,  with  a  response  capable  of  being  analysed 
by  the  methods  of  the  theorv  v.f  elastic  wave  propagation  and,  .Hi  frac¬ 
tion. 

A  transducer  to  measure  ground  stresses  should  have  the  following 
char  ante r is t  ics : 

1,  The  strength,  to  wit  hat  and  high  stresses 

2 .  Elastic  properties  to  produce  measurable  strains 
1.  Fve  of  instal  lat  Lor.  and  recove  rab  1 1  l  tv 

d.  Insensitivity  to  orientation 

f.  Small  rise  time,  minimized  overshoot  ,  and  large  internal  tran¬ 
sient  den  ,lv 

»••.  Compatibility  with,  standard  aux ilia,  c  equipment 


Pie  transducer  w 

e  propose  t 

examine  is  a 

r.  clast 

ie  spherical  in¬ 

elusion  to  be  buried 

in  the  gtour.d 

!,  and  sens  it 

1  v  <.■  ,  fh 

rough  strain  gage* 

to  *-he  mean  •■tress  or 

pressure  at 

its  origin. 

Such  a 

device  is  ins,  r.- 

sitive  to  its  angular 

or  lent  at  ion 

hot:-,  because 

Of  it  3 

svmnetrv  and  be- 

cause  of  the  scalar  n 

atvire  of  the 

mean  stress 

or  pressure.  As  our  fir 

task  we  will  examine 

the  response 

chara . t er  i  s  t 

i  o  s  o  f 

.such  an.  i  r. c  1  us  i  on 

task  we 


n  meet  tie  requirement  of  high  strength,  the  material  of  the  sphere- 
will  most  likely  have  elastic  properties  different  from  those  of  the 
ground  material.  Hence,  diffraction  effects  are  likely,  and  their  ef- 
fec*.  on  the  response  of  the  transducer  should  be  examined.  Toward  this 
end  we  treat  the  problem  of  the  diffraction  of  a  transient  compressions! 
pulse  by  an  elastic  isotropic  spherical  inclusion  bonded  to  an  infinite 
elastic  isotropic  matrix.  Then  some  design  criteria  are  presented  for 
the  selection  of  the  material,  size,  ■'d  response  times  of  a  spherical 
strain  transducer. 
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II.  KPE 


RESPONSE  GF  AN  ELASTIC  INCLUSION 


The  study  of  the  response  of  an  elastic  inclusion  has  received 
much  attention  in  the  last  decade.  These  studies  have  dealt  mainly 
with  harmonic  compressions!  waves.  ^  ^  Recently,  however,  Mow^’^ 
has  examined  the  transient  behavior  of  a  rigid  inclusion  and  has  cor¬ 
related  the  separate  displacements  of  the  inclusion  and  the  ground. 

He  found  that  ’he  inclusion  and  free-field  displacements  could  diverge 
at  early  times  even  for  equal  densities,  satisfying  an  essential  cri¬ 
terion  of  field  testing.  Similarly  the  stresses  in  an  elastic  inclu¬ 
sion  are  expected  to  differ  in  both  magnitude  and  phase  from  those  in 
the  free  field.  The  extent  of  this  divergence  and  the  factors  tht 
will  minimize  it  are  discussed  below. 


BASIC  EQUATIONS 

(3) 

Since  the  harmonic  solution  has  been  discussed  elsewhere,  we 
will  only  briefly  review  the  relevant  equations. 

We  suppose  both  the  Inclusion  and  matrix  to  be  homogeneous,  Iso¬ 
tropic,  linearly  elastic  mediums  in  which  the  usual  stress-strain  re¬ 
lationships  held.  For  each  material  there  are  two  elastic  constants 
and  two  wave  speeds:  a  compresslonal  or  longitudinal  wave  c^ ,  and  a 

shear  or  transverse  wave  c.  . 

s 

Outside  the  inclusion  we  Imagine  a  source  of  either  plane  or  spher¬ 
ical  harmonic  compresslonal  waves  with  a  displacement  potential  given  by 


p  iuifz  cdjt)/cdl 

*  =  V 


.s  V 


iu)(R-cdlt)/cdl 


plane  wave 


spherical  wave 


u  " 


where  z  is  an  axial  coordinate  (Fig.  1),  R  is  the  distance  from  the 
source  to  the  field  point,  and  u  is  the  displacement.  The  potentials 


r’/v. . 


Geometry  of  spherical  inclusion 
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represented  by  Eqs.  (1)  and  (2)  give  the  displacement  in  a  medium  with¬ 
out  an  inclusion.  The  presence  of  the  inclusion  will  scatter  waves  into 
the  medium  and  refract  waves  into  the  inclusion.  Tne  scattered  and  re¬ 
fracted  fields  will  contain  both  compressions!  and  transverse  waves. 
Inside  the  inclusion  the  motion  is  of  the  form  of  standing  waves 

or  vibrations,  and  the  refracted  displacement  potentials  take  rhe 

.  (3) 

form, 


u  =  V4>  +  V  x  fe^Y/ao) 


*  “  -£  Cnjn(u)r/cd2)  Pn  (cos  0) 
n=0 


Y  =  - y,  D  j  (wr/c  „)  P  i COS  0) 
*-■'  n  n  s2  n 
n^O 


(4) 


where  j  (a)  are  spherical  Bessel  functions  and  P  (x)  are  Legendre  poly- 

n(h)  n 

nomials . 


AVERAGE  AND  CENTRO IDAL  STRESSES  AND  DISPLACEMENT 

The  induced  stresses  are  derived  from  Eq .  (4)  using  the  generalized 

(3) 

Hooke's  Law  and  have  been  tabulated  by  Pao  and  Mow.  The  stresses  and 
displacements  depend  on  all  the  excited  modes  in  the  inclusion.  However, 
if  the  inclusion  is  being  used  as  a  transducer,  then  knowledge  of  the 
stresses  at  one  point  will  be  sufficient.  The  most  logical  point  to  ex¬ 
amine  is  the  origin. 


Using  the  smal 1 -argument  limit  for  the  spherical  Bessel  functions, 


(6) 


Vu)aK)  1.  3.5  .  .  .  (2n  +  1) 


the  d  placement  and  stresses  at  the  origin  become 


u 

z 


3c 


fn. 


dl 


2W 


(5) 


T  -  -P  +  2S 
zz 


1  ■  T  *  -P  -  S 

xx  yy 


where  k,,  ■  t^e  rat*°  compresaional-wave  speed  to  shear 

wave  speed  in  the  inclusion.  Thus  the  displacement  at  the  origin  de¬ 
pends  only  on  the  first  dilatatlonal  and  shear  modes.  The  mean  pres¬ 
sure  P  at  the  origin  depends  only  on  the  symmetrical  modes,  and  the 
maximum  shear  stress  3S/2  depends  solely  on  the  vibratory  modes  asso¬ 
ciated  with  C2  and  D^. 

Similar  remarks  apply  dor  the  ave*.  'ge  displacement  or  motion  of 
the  center  of  mass  of  the  inclusion  and  average  pressure  over  the  en¬ 
tire  sphere,  i.e.  , 


<u>  • 


-  2DlJ]L(62) 


(7) 


<P>  ~  <~B2V*u> 


-B, 


"d2 


2j1(a2) 


(8) 


where  a„  ■  ca  "  cua/c._  and  0„  -  wa/c 

2  dZ  2  s2 

It  is  evident  that  the  behavior  of  either  the  stresses  at  the  or¬ 
igin  or  the  average  mean  pressure  depends  only  on  the  frequency  response 
of  the  coefficients  Cq,  C2,  and  D2<  These  coefficients  are  determined 


_7.* 


from  Che  boundary  conditions.  For  a  bonded  inclusion,  the  boundary 
conditions  require  continuity  of  stress  and  displacement,  i.e., 

u(s)  +  u(I>  -  u(r)  ] 
r  r  r 


(s)  ,  (I) 

u  +  u„ 

6  u 


on  r  *  a 


T<S)  +  T(I) 

rr  rr 


T(s)  +  t(D  _  T(D 

r6  rO  re 


where  s,  I,  and  r  indicate  "scattered,"  "incident,"  and  "refracted," 

* 

respectively.  For  our  purposes,  it  is  sufficient  to  know  that  C0  de~ 
pends  only  on  the  first  and  third  boundary  conditions  of  F.q.  (9).  That 
is,  the  pressure  at  the  origin  or  the  average  pressure  ever  the  entire 
inclusion  does  not  depend  on  whether  the  tangential  displacements  or 
shear  stresses  at  the  Interface  are  continuous.  Thus  the  same  value 
of  Cq  will  result  if  "slipping"  boundary  conditions  of  zero  shear  stress 


on  the  interface  are  applied. 


(*)  .  (T) 

u  +  u 
r  r 


T(s)  +  t(D 
rr  rr 


■  0 


on  r  *  a 


0  -  T 


The  induced  pressure  at  the  center  of  our  transducer  inclusion  will 
thus  be  insensitive  to  the  bonding  or  grouting  as  long  as  the  rratrix 


These  equations  are  given  in  a  more  explicit  forrr.  in  Ref.  1 . 
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m&kes  normal  contact  with  the  sphere.  For  this  reason,  we  study  only 
the  frequency  and  transient  response  of  the  mean  pressure  at  the  ori¬ 
gin  of  the  inclusion.  The  coefficient  Cq  is  given  by  the  expression 


where 


co  "  -N0iaeia<|  ^  a<a  + 


[“kr1  * . •■] 

i,(ca)j 

plane  wave 

(11) 


No-  < 


50  H  koH) 


spherical  wave 


where  1  is  the  distance  from  the  wave  source  to  the  inclusion  center, 
and  h^(x)  is  ti  -  zero-order  spherical  Hankel  function. 

Let  us  next  examine  the  consequences  of  choosing  either  a  plane 
or  a  curved  incident  compressional  wave.  The  mean  pressure  due  to  the 
incident  weve  is  given  in  the  two  cases  by 


F  -  -B  V*uvI)  - 
c  1  -  1 


_®1 

c2  at2 

ax 


(12) 


where  $  satisfies  the  wave  equation.  ^  Substituting  Eqs.  (1)  and  (2) 
into  the  above  gives 


plane  wave 


(13a) 


’  “  Bi  (JL-\  3 

C  1  \Cdl  /  I 


*5  h0  (*) 


where  hQ(x)  •  e*X/ix. 


spherical  wave 


(13b) 
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Comparlng  Eqs.  (13a)  and  (13b)  ,  ^  we  find  that  the  relationship 
between  the  incident  pressure  and  inclusion  pressure  is  the  same,  i.e., 


,  n  c  iae 
P(a)  =  b 


-la 


c  B  H(a) 


A(a)P 


(14) 


H(a) 


^  a (a  +  i)j0(ca)  - 


u  +  i)  +  ia2 


j  x  (ca)/ 


We  need  no  longer  distinguish  between  spherical  or  plane  waves,  since 
the  response  is  identical  in  both  cases.  The  same  can  be  said  for  the 
displacement  and  shear  stress  3.  In  fact,  as  long  as  a  quantity  de¬ 
pends  on  only  one  mode,  the  relation  between  its  value  inside  and  its 
value  in  the  free  field  remains  the  same  whether  the  incident  wave  is 
plane,  spherical,  or  generally  dilatational .  However,  for  a  quantity 
such  as  t  ,  the  axial  stress  which  depends  on  the  zero^  and  second 
modes,  this  is  not  true;  i.e.,  if 

i  (r  =  0)  =  t*  g  (u;)  plane  wave 

E  Z  —  Z  Z  p 

t  (r  -  0)  »  g  (oj)  spherical  wave 

zz  -  zz  s 


then  8p(w)  i  g^(u).  This  gives  further  Impetus  to  studying  the  mean 
pressure  P  rather  Lhan,  for  example,  the  axial  stress  t 

z  z 

Examining  the  frequency  response  for  the  mean  pressure  at  the  or¬ 
igin,  we  observe  that  at  low  frequencies,  i.e.,  -*•  0,  n  *■  0 

J0(cn)  -  1  ,  jj(ci)  ‘  ,.»/3 


P 


c 

B 


1 


(lo) 


This  limit  yields  the  same  pressure  as  the  static  problem  of  a  spheri¬ 
cal  elastic  inclusion  embedded  in  an  elastic  matrix  with  hydrostatic 


-10- 


pressure  P  at  infinity. '7^  The  ratio  of  output  to  input  P/P  is  a 
c  c 

constant  at  low  frequencies  and  can  be  called  an  amplification  or  gain 
factor  Aq.  The  gain  factor  is  presented  in  Table  1  for  different  in¬ 
clusion  materials  in  a  granite  matrix.  At  high  frequencies 

P(a)  -*■  P  e*a(pc  cos  ca  -  i  sin  ca)  ~  (16) 

a—  C  B 

The  product  pc  is  called  the  impedance  ratio,  and  its  departure 
irotn  unity  gives  a  clue  to  the  nature  of  the  reflection  of  dilatation 
waves  from  the  boundary  separating  two  media.  If  pc  -  1,  the  gain  is 
again  constant,  whereas  if  the  impedance  ratio  departs  from  unity,  the 
denominator  has  periodic  maxima  and  minima.  Thus  the  gain  exhibits 
resonance  peaks  with  frequency,  or  the  internal  reflections  reinforce 
one  another  at  certain  frequencies.  The  resonances  at  large  frequency 
are  easily  found:  for  pc  >  1,  cn  ■  ns/2,  and  n  an  odd  number, 

(P/P  )tnax  -  p(c)^/B  (17a) 

c 

and  for  pc  <  1,  ca  •  nir/2,  and  n  an  even  number, 

(P/P  )ma:  -  c2/B  (17b) 

c 

The  frequency  response  curves  using  the  exact  relation  Eq .  (14) 

* 

are  shown  in  Figs.  2  through  5.  They  show  the  resonances  at  the  proper 
frequencies  when  a  is  large.  Values  obtained  previously  check  the  com¬ 
puted  maxima  of  these  curves. 

An  ideal  transducer  is  one  for  which  the  gain  or  amplification 
factor  is  independent  of  frequency.  Thus  for  measurement  of  elastic 
waves  in  granite,  aluminum  or  titanium  would  stake  a  better  transducer 
then  plastic  or  lead,  for  example  (Figs.  2  through  5).  Furthermore, 
we  can  cay  that  unleae  the  matrix  and  inclusion  transducer  are  ’'matched" 
CjSf  -  1),  the  internal  reflections  can  mask  the  incident  signal. 


Por  discussion  of  Cases  A  end  B  in  Fig.  2,  see  page  20  and  Table  2. 


■8 
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f 


High-frequency  limit  1/AqB 


Fig. 2—  Frequency  response  of  the  pressure  at  the  center  of  a  sphericol  elastic  inclusion  in  an  infinite 

matrix;  effect  of  mismatch  of  shear  modulus 


Wave  number  a  =  ao/c^ 

Frequency  response  of  the  pressure  in  a  lead  inclusion  in  a  granite  matri; 


Fig. 4 —  Frequency  response  of  the  pressure  in  <j  polyethy lene  inclusion 

in  a  granite  matrix 


Wave  number  a  wa  e, 

dl 


Frequency  response  of  the  pressure  in  a  titanium  inclusion 
in  a  granite  matrix 
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111.  TRANSIENT  RESPONSE  OF  AN  EMBEDDED  ELASTIC  INCLUSION 

DIRECT  PROBLEM 

If  the  incident  signal  or  pressure  is  nonperiodic,  the  inclusion 

pressure  at  the  origin  may  be  obtained  a3  an  inverse  Fourier  transfo^... 

of  the  product  of  the  gain  A(a)  (Eq.  (14))  and  the  transform  of  the  in- 
(9) 

cident  pressure,  i.e,, 


P(t) 


A(a)F[P  (t) ]e_laTda 
c 


> 


F[P  (t) ]  =  f  P  (T)eiaTdi 
c  J  c 


(18) 


where  P  (x)  would  be  the  pressure  at  the  origin  in  the  absence  of  the 
inclusion.  The  time  t  is  normalized  by  the  time  the  incident  wave 
transits  half  the  inclusion,  i.e.,  a/c^,  and  a  is  a  nondimens ional 
frequency. 

Suppose  P  (t)  has  the  form 


pc(x) 


r 

o, 

PQf(T), 


T  <  0 


T  >  0 


A  x 

where  f(r)  -*•  e  is  t  +  ®.  In  units  of  the  half  transit  time  a/c^j, 
the  pulse  arrives  at  the  inclusion  surface  at  x  ■  -  !•  time  for  a 

compressional  wave  to  transit  the  distance  "a"  in  the  inclusion  mate¬ 
rial  is  c.  We  should  then  expect  the  response 


P (x)  -  0 


1  r  i  <  -1  +  c. 


Tc,  show  that  this  is  the  case,  the  calculus  of  *-esiduee  is  used 
to  evaluate  Eq .  (18)  and  to  show  that  all  poles  rf  the  integrand  lie 
in  the  lower  half  of  the  complex  plane.  (See  Fig.  *».)  In  the  exami¬ 
nation  of  the  function  A(a)  in  the  com, lex  a-plane,  we  note  that  the 
singularities  are  determined  by: 


Real  H(a)  »  0 


Im  H(a)  *  0 


(19) 


It  can  be  shown  that  the  roots  of  Eq.  (19)  always  come  in  pairs, 

*  t 

i.e.,  if  a  is  a  root,  -a  is  also  a  root.  To  show  that  the  imaginary 
part  of  the  root  is  always  negative,  we  note  that  as  real  a  where 

a  »  ft  -  iy,  the  equations  for  the  roots  of  Eq.  (19)  take  this  form 
(where  H(a)  is  given  by  Eq.  (14)): 


cos  cft(pc  cosh  cy  -  sinh  cy)  ”  0 


sin  cft(cosh  cy  -  pc  sinh  cy) 


ThuJ  the  roots  are  determined  by 


Case  1:  pc  <  1,  sin  c!l  *  0,  eft  ■  nir/2,  n  even 


tanh  cyQ  -  pc,  y  -  yQ 


(20a) 


Ce  2:  pc  >  1,  cos  Cft  *  0,  eft  *  nn2,  n  odd 

(20b> 

tanh  cYq  “  1/pc 

For  high  frequency,  the  radiation  damping  coefficient  y  approaches  a 
positive  constant  if  JTF  4  1. 

+ Asterisk  indicates  conjugate. 
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Since  I-  a  <  0,  there  are  no  poles  in  the  upper  half  plane,  and 
substitution  of  Eq.  (14;  for  A(a)  into  Eq.  (18)  reveals  that  P(t)  «  0 
for  t  +  1  -  c  <  0,  which  agrees  with  our  initial  supposition. 

For  Im  a  <  0 


A(a)  n,  e 


i(l-c)a 


and  P(t;  for  (t  +  1  -  c  >  0)  is  given  by  the  sum  of  the  residue  of 
A(a)F[f]  in  the  lower  half  plane. 

For  the  special  case  of  a  step  function, 


f 


->  T 
e 


\  -*■  0 


it  follows  that 


P(t)  -  { 


0  ,  T  <  -1  +  C 


T  -(y  +i-Q  )  (t+1  ) 
A0  -  1  - 


(21) 


T  >  - 1  +  C 


where  the  sum  is  over  the  conjugate  roots  of  H(a  )  ■  0  and  the  property 

H'(-a  )  »  -H  (a)  has  been  f’ken  into  consideration.  For  large  j  the 

roots  take  the  form  given  in  Eq.  (20). 

The  first  term  in  Eq .  (21)  is,  to  within  a  factor,  the  incident 

pressure.  The  second  is  a  transient  pressure  due  to  diffraction  and 

internal  reflections .  The  maximum  rate  of  decay  of  the  transient  is 
-y  T 

proportional  to  e  0  where  y^  is  given  by  either  Eq .  (20a)  or  Eq .  (20b). 
The  normalized  "damping"  time  is  l/y^.  For  a  granite  matrix  and  alumi¬ 
num  inclusion  this  time  is  1.5  half  transit  times,  whereas  for  a  plastic 
inclusion,  l/y^  “  30  half  transit  times. 

The  complete  solution  may  be  obtained  by  using  the  first  few  roots 
of  H(a)  **  0,  found  numerically  until  y ^  -*•  y^.  The  remaining  terms  in 


t 


Prime  sign  indicates  derivative. 


-20- 


the  series,  using  the  approximation  Eq.  (20)  for  a  ^ ,  constitute  a  Fou¬ 
rier  series  and  may  then  be  summed  analytically. 

As  an  illustration,  consider  the  cases  for  which  the  frequency  re¬ 
sponse.  is  calculated  in  Fig.  2  (the  upper  curve  is  close  to  that  for 
an  aluminum  inclusion  in  a  granite  matrix) .  In  these  cases  the  densi¬ 
ties  and  compressional  velocities  are  matched,  but  the  shear  speeds  of 
inclusion  and  matrix  are  different,  i.e.,  pel  -  1,  p  f  1.  The  first  few 
roots  and  residues  for  each  case  are  presented  in  Table  2,  We  observe 
that  the  magnitudes  of  th ’  third  and  fourth  modes  are  about  one-half 
of  the  first  symmetric  vibration  mode  of  the  inclusion  in  the  matrix. 
The  relative  importance  of  the  higher  modes  was  also  found  by  Skalak 
and  Friedman in  studying  the  reflection  of  an  acoustic  step  wave 
from  an  elastic  cylindrical  shell.  We  note  also  that  the  real  parts 
of  the  roots  came  close  to  the  peaks  in  the  frequency  response  curves 
(Fig.  2). 

The  pressure  response  to  an  incident  compressional  step  wave  at 
the  center  of  an  aluminum  inclusion  in  a  granite  matrix  is  given  by 

P(t,  -  PQ(22/9)[l  +  (1.301)e_1'251(T+I)sin[3.532(T  +  1)  +  9, 

-t-  '0.8465)e_1'54°(T+1)sin|6.839(T  *  1)  +  02j 
+  (0.6827)e-1‘ /19(T41)sinfl0.04(';  +  1)  +  83j  +  ...j 

for  t  >  0. 

This  curve  is  plotted  in  Fig.  7,  using  only  the  first  three  modes. 
(Note  that  since  pc  •  1,  we  cannot  use  Eq .  (20)  to  analytically  sum  the 
rest  of  the  terras  in  the  series.)  Even  though  the  series  converges 
slowly  tor  any  given  time,  the  exponential  radiation  damping  terms  en¬ 
sure  that  the  transient  will  be  negligible  in  a  few  transit  times  for 


this  case. 
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LARGE-TIME  AND  SMALL-TIME  RESPONS E 

Let  us  now  examine  the  behavior  of  Eq.  (18)  for  x  +  ®  or  its 
equivalent  a  -  0.  For  this  case  A(a)  -*•  A(0)  5  AQ,  a  constant,  and  the 
remaining  integral  is  the  inverse  of  F[f] 

P ( t )  AnP  (t  +  1)  (22) 

U  c 

T-MC 


(It  should  be  noted  that  the  phase  of  the  right-hand  term  is  arbitrary 
for  t  ®. ) 

For  small  times  we  examine  Eq.  .14)  under  the  limit  a  ■*  for 
which  Eq .  (18)  becomes 


P  (t) 


,->.3  i  r°°  -ta(T+l)  .  , 

-  (c)  1  /  e _  F[  f  ]  da 

^  o  2tt  J  (JTC  cos  ca  -  i  sin  ca) 

u  —  Oc 


(23) 


This  may  be  inverted  if  we  rewrite  Eq.  (23)  in  the  form  of  a  Laplace 
transform,  setting  a  *■  is 


3  .1-  s(T+l)L[f]ds 

P(T)  *  PAP  —^7  / 

~  A 


'0" 


B 


,  (Pf  cosh  cs  +  3inh  cs) 

-  i  a> 


(24) 


or 


—  3  ■  i 

P(t)  -  Pnd(~'|  sinh  Zy  ™  /  — 

0  \  b  I  Q  2-i  J  .  co 


i-  • ( T+l) 


L [ f ) ds 


.  cosh  c(s  +  yn ) 
i »  0 


(24a) 


where  t  nh  cyq  «  i/pc',  pc  >  1  and 


3  r  i-  s(t+1) 

P(t)  -  P  o  cosh  cy  /  -2? - VL/ JL*3.* 

K  ’  r0o  _  cosn  cy0  2wl  ^  einh  g(f  +  ) 


( 24b) 


when  tanh  cYq  “  pc,  pc  <  1.  The  integrals  in  Eq* .  (24a)  and  (24b)  can  be 


2  r 


evaluated  by  expanding  ej.:her  1/coah  i  or  1/ainh  r  in  powers  of  e  , 
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from  which  we  see  that  Eq.  (24)  is  a  series  of  modulated  step  functions 
separated  in  time  by  2cn.  These  represent  reflections  inside  the  in¬ 
clusion  . 

Integrals  of  the  form  Eq.  (24)  may  be  found  in  a  table  of  inverse 
Laplace  trail.-  forms .  ^  Finally  we  have,  for  small-time  response, 

Case  1  ("soft"  inclusion):  pc  ■■  1 


P(t) 


2  sinh  cy 
B 


0 


L(-Dn 

n-0 


- (2n+l) cyq 
e 


f[:  +  1  -  (2n  +  l)c]u[t  +  1  -  (2n  +  1 ) c ]  (25a) 

Case  2  ("hard"  inclusion) :  (Tc  <  1 


P(t) 


cosh 


CY  , 


3'  -(2n+l)cy 

v >  0 

L  e 

n*-0 


f[t  +  -  -  (2n  -*■  lk]u(:  +  1  -  (2n  +  l)c]  (25b) 

where  u(t)  is  the  unit  step  function. 

These  expressions  are  applied  to  a  few  examples  as  shown  in  Figs. 
8,  9,  and  10  for  the  case  of  an  Incident  step  pressure.  For  Case  1 
(TFcT  y  1)  the  inclusion  acts  as  an  underdamped"  transducer  (Figs.  9  and 
10)  and  as  an  "overdamped"  transducer  for  Case  2  C'c  <  1)  (Fig.  q/. 

The  damping  is  due,  of  course,  to  the  fact  that  part  of  the  reflected 
waves  is  radiated  out  of  the  inclusion  into  the  matrix. 

If  we  t  rv  to  extend  Eq .  (25)  for  large  time  for  the  case  of  the 
step  Input.,  we  obtain  in  both  cases 

Pit)  -  5"  P,  (26) 

R 


where  we  have  used  the  Met  that 


transient  response  of  lead  inclusion  in  a  granite  matrix 
step  wave;  decay  of  internal  reflections 


granite  matrix  for  an  incident  compressiona!  step  wave;  decay  of  internal  reflections 
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10  -  l'2n+i)  cy  _ 

2l«  '  ° 

n«0 


sinh  cy„’ 
u 


“  „  ~(2n+i)cy 

2E  (-!)”«  ° 


1 


n»0 


cosh  cy 


0 


Thus,  1  ;r  the  step  input  the  short-time,  asymptote  given  by  Eq .  (26) 
using  the  approximation  Eq.  (23)  differs  from  the  long-time  approxima¬ 
tion  Eq.  (22;  by  a  constant.  This  difference  will  be  discussed  shortly. 
The  above  factor  represents  the  asymptotic  maximum  (minimum)  in  the  fre¬ 
quency  response  for  the  case  pc  <  1  (pc  >  1). 

For  the  case  Jc  >  1 


c2  /  \  c2  ~2Zy0 

!'  =  —  (  2  cosh  cy  e  -  1  )  =  — •  e 

B  '  0  •  B 

is  a  measure  of  the  amount  of  overshoot  for  early  time. 

We  see  that  this  formulation  has  the  advantage  over  the  Fourier 
representation  Eq,  (18),  since  for  a  given  time  Eq .  (25)  contains  a 
finite  number  of  terms.  In  the  latter  expression  1/C y^  again  plays 
the  role  of  a  "transient  damping"  time. 

For  the  case  of  equal  impedances  pc'  *>  1,  we  return  to  the  expres¬ 
sion  Eq.  (13).  The  small-time  expression  is  given  simply  by 

_2 

P(t)  *  •=-  Prf(x  +  1  -  c) 

B  0 

The  small-time  solution  Is  exactly  analogous  to  the  penetration  of 
a  compresslonal  wave  into  an  infinite  acoustic  barrier  of  material  dif¬ 
ferent  from  that  of  the  matrix  and  normal  to  the  incoming  wave.  For 
that  problem,  if  the  acoustic  impedance  ratio  is  unity,  there  are  no 
reflections,  but  rather  a  delay  in  the  signal,  which  is  what  we  found 
previously. 
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IV.  THF,  SOLUTION  OF  THE  INVERSE  PROBLEM 


While  the  goal  of  the  analyst  is  to  predict  the  output  or  response 
of  a  transducer  given  the  input,  the  object  of  the  experimentalist  is 

the  reverse.  That  is,  the  latter  wishes  to  solve  the  inverse  problem 

* 

of  finding  the  input  given  the  output  of  the  transducer. 

THE  METHOD  OF  FOURIER  TRANSFORM 

Returning  to  Eq .  (14),  we  write  an  expression  for  the  free-field 
pressure  in  terms  of  the  mean  pressure  P  at  the  origin  of  the  in¬ 
clusion 


P  (T) 

c  c  2it 


e-icx(T-l)Fjpj 


X 


(-i)(a  +  i)j0('ca)  - 


4(1  -  y) 
2 

ykl 


(a  +  i)  +  ict 


(27) 


One  may  readily  perform  this  inversion  if  the  integral  form  of 
Jn(e)^  is  used,  i.e.  , 

VEa)  “  i  /x  elcatdt 

J2(ea)  -  -  ^  J  ei^at(3t2  -  l)dt 


and  if 


3J1(r)  -  zj^J0(e)  +  J2(*)j 

The  inversion  is  performed  by  switching  the  order  of  integration 
and  making  the  necessary  assumptions  of  integrability .  Thus 


|  See  Ref.  11  for  an  example  from  quantum  mechanics. 

! 

i 
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l'c  (x)  »  ~  f  A 

c  2c  J  .  pc 


dP 

<3t 


(t  -  1  -  ct>  +  P(t  -  1  -  St)  dt 


_  B  f  J  4  (1  -  ^ 


uK 


dP 

dr 


(t  -  1  -  ct) 


p(t  -  1  -  ct)l  -  ~  (t  -  1  -  ct»-(l 

J  dr" 


t2)dt 


Integrating  by  parts  and  collecting  terms,  the  above  expression  assumes 
the  form 


P  (t)  -  — L  |(1  +  A  P(t  -  1  +  c)  +  (l  -  ~)p(t  -1-5) 
C  2(E)  L  Pc  \  TKJ 


)  -  f)  z  L  p(t  -  1  s)ds 


+  ^—"4  [  (s2  +  2s  -  c2)P(t  -  1  -  s) ds  (28) 

(?)  -c 


Note  that  ifu*c*8«p«l,  we  obtain  the  identity  Pq(t)  "  P(t). 
For  long  time  or  t  -  1  »  2c,  we  may  remove  P(t  -  1  -  s)  from  the  in¬ 
tegrands  and 


PC(T) 


I  LI  -  A  Cl  -  p> 

5  Lpc  3  Pk2 


P(t  -  1) 


Aq1P(t  -  1) 


(29) 


which  is  vhat  we  obtained  for  the  direct  problem,  Eq.  (22). 

It  is  interesting  to  note  that  Mow^  found  an  exact  expression 
for  the  inverse  problem  of  a  rigid  sphere.  However  in  that  problem  the 
incident  displacement  was  found  in  terms  of  derivatives  of  U  as  well 
as  an  integral.  Equation  (28)  involves  only  weighted  integrals  for  the 
incident  stress.  This  operation  could  be  built  into  an  electronic  "black 
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box."  The  recorded  signal  P(t),  considered  as  an  Input,  could  be  de¬ 
layed  and  Integrated  electronically  to  produce  an  output  proportional 
to  the  original  incident  pressure.  The  operations  in  Eq.  (28)  thus  de¬ 
code  the  transducer  output. 

The  solution  of  the  inverse  problem  for  the  center  pressure  was 
possible  because  CqA(u)  in  Eq.  (6)  was  a  linear  combination  of  spheri¬ 
cal  Bessel  functions.  No  such  simple  relation  results  for  the  s..ear 
stress  S,  Eq.  (6),  or  the  displacement  of  the  center  of  mass  <£>,  Eq. 
(7). 


THE  METHOD  OF  SUCCESSIVE  REFLECTIONS 

Having  obtained  an  exact  expression  for  the  incident  pressure  P  (t) 
in  terms  of  the  center  pressure  in  the  inclusion,  we  can  then  look  at 
Eq.  (28)  as  an  integral  equation  for  P(t)  in  terms  of  the  incident  pres¬ 
sure.  We  will  in  fact  chow  that  by  this  method  an  exact  solution  for 
P(t)  can  be  obtained  by  a  finite  time,  step-by-step  algorithm. 

To  l.0in,  we  shift  the  time  by 


Tj  =  t  -  1  +  c  and  n  =  -  c  -  s 

so  that  Eq.  (28)  takes  the  form 


P  (t  +  1  -  l)  -  ~~~2 
2(c) 


i1  *  +  (!  -  -  2« 


-  (i  -  AH  /  pc 

'  p  '  c  J  _ 


n)dn 


+  _L_  fr  i 

3  2 

(c)  Okj^  Jt1-2C 


P  (n ) 


(t. 


n  -  c) 


2 


+  2 (t^  -  n  -  c) 


dn 
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We  assume  that  the  center  pressure  is  zero  until  *  0,  or  that 
the  Incident  pulse  arrives  at  the  origin  at  time  -  1  -  c.  These 
conditions  imply  the  following: 

P^)  -  0  Tj  <  0 

P  (t.)  -  0  t,  <  1  -  c 

Cl  -l 

P(t^  -  2c)  *  0  t ^  <  2c 

,0 

I  P(n)dn  -  0 

v25 

Since  c  is  the  half  transit  time  in  the  inclusion,  the  center  of 
the  inclusion  does  not  see  an  internal  reflection  until  time  -  2c. 
Thus  for  0  j<  <_  2c, 


F 

c 


U, 


i-  x 


c) 


[E(1  +  sf)p<V 


n)P(n)dn 


(30) 


where 


K(t1  -  n) 


+  2(t1  -  n 


This  is  an  integral  equation  of  the  Volterra  kind.  Since  K(x^  -  n) 
is  a  polynomial,  Eq.  (30)  may  be  solved  by  successively  differentiating 
it  with  respect  to  t^,  using  the  identity 


h(t1  -  n)P(n)dn 


K(0)P(t1)  + 


n)P (n)dn 


This  results  in  an  ordinary  linear  differential  equation  with  con¬ 
stant  coefficients  for  P(t^),  with  initial  conditions  determined  at  each 
differentiation  of  Eq.  (30),  i.e.. 
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P  (t.  +  1  -  c) 
c  1 


2(c)3  L 


;(l  +  ==)p'(t^  +  K(0  )P(T1) 


K#(0)P(t1)  +  K'(0)P<t1)J 


with 


(31) 


P  (1  -  c)  -  ~ 2-T 

C  2(c)2 


P  (1  .  e)  -  -2-j 

/rr\  J 


Pc(l  -  c) 


2(c) 


2(c) 


(l  +  #(ot> 

c(l  +  =)p(0+)  +  X(0)P(0+)j 


[E(1+#(0+ 


)  +  K(0)P(0+)  +  K/(0)P(0+) 


Ola) 


(31b) 


Ole) 


For  our  problem 


K(°)  .(i-  l) 

b kj 


(0)  -  (1  _  c) 

bk^ 


/  (o)  -  ^ 

yki 


The  solution  consists  of  a  particular  solution  determined  by 

P  (t,  +  1  -  c)  and  a  transient  solution  of  the  form 
c  1 


3  ib  t 

P(V  -  E  J  0  <  X.  -  2c 

J-]  J 


(32) 


where  the  roots  of  the  following  algebraic  equation  are  represented 
by  ty 


I 
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c(l  +  =)(18)3  +  K(0)  (if?)2  +  K'(0)(it?)  +  K'(0)  «  0  (33) 


The  three  constants  a^  are  determined  by  the  initial  conditions,  Eqs. 
(31a),  (31b),  and  (31c). 

Once  the  solution  is  known  in  the  interval  0  <  £  2‘c,  the  solu¬ 
tion  may  be  extended  into  the  interval  2c  <  <_  4c,  since  P(t^  -  2c) 

will  be  known,  at  will 


r 


V2c 


P (n) K(t1  -  n)dn 


A  comparison  between  the  solutions  obtained  in  Eq.  (21)  (method 
of  residues)  and  Eq.  (32)  (method  presented  here)  reveals  a  different 
in  the  number  and  magnitude  of  eigenvalues  present  in  each.  In  Eq .  (32), 
there  exists  a  finite  number  8^;  in  Eq.  (21),  there  is  an  infinite  num¬ 
ber  Oj.  This  can  be  explained  by  the  fact  that  the  two  forms  of  the 
solution  are  valid  over  different  time  domains.  Equation  (21)  i3  valid 
for  all  time  and  must  contain  information  about  the  periodic  discontin¬ 
uities  in  pressure  due  to  reflections.  The  method  discussed  here  ob¬ 
tains  a  solution,  Eq.  (32),  valid  only  between  reflections  (2nc  <  < 

2(n  +  l)c)  and  continuous  in  these  intervals  as  long  as  P  is  continuous. 

We  may  note  that  this  method  may  be  applied  to  the  conversion  of 
certain  integrals  of  the  form 


f(t) 


e  ^at~H(a)d3 

E  «V«> 

n— M 


where 


H( 


c  °° 

a)  ■  J  e*ath(t)dt 


and  the  g  (a)  are  periodic  functions  and  bounded.  One  must  also  assume 
N  n 

that  a  g^a)  ”  0  has  no  roots  on  the  real  axis, 

n* -M 


V.  SUMMARY  OF  RESPONSE  CHARACTERISTICS 


We  now  recapitulate  the  main  results  of  the  various  methods  used 
to  analyze  the  response  of  an  embedded  elastic  spherical  transducer. 

1.  The  pressure  or  mean  stress  response  at  the  center  of  the  in¬ 
fusion  will  be  independent  of  the  orientation  of  the  inclusion  relative 
to  the  wave  front  and  will  only  depend  on  the  spherically  symnetric  modes, 
as  shown  in  Eq .  (6). 

2.  The  pressure  response  at  the  center  of  the  inclusion  will  not 

depend  on  whether  the  incident  wave  is  plane  or  spherical.  The  stress 

7  ,  however,  will  be  sensitive  to  the  spherical  character  of  the  wave 

zz 

front . 

3.  The  higher  free  spherical  vibrational  modes  of  an  embedded  in¬ 
clusion  are  characterized  by  coflanensurable  frequencies  and  radiation 
damping  independent  of  frequency. 

A.  The  rate  of  decay  of  internal  reflections  will  depend  on  the 
impedance  ratio  pc'  and  will  have  a  characteristic  time  y  ^  given  by  Eqs. 
(20a)  or  (20b)  (see  Table  1). 

5.  The  pressure  response  for  su-ill  time  is  characterized  by  behav¬ 
ior  which  ia  either  "underdamped"  (Eq .  (25a))  or  "overdamped."  For  the 

iorraer,  the  overshoot  above  the  small-time  asymptote  is  given  by 
2  2  c  y  ~ 

(c)  e  'O/B  for  an  Incident  step  wave  (see  Table  l/. 

6.  The  inverse  problem  may  be  solved  for  the  incident  pressure  in 
terms  of  the  pressure  st  the  certer  of  an  embedded  inclusion,  from  which 
an  exact  solution  may  be  obtained  for  the  direct  problem  by  a  reflection- 
by-reflection  method.  His  is  not  possibls  for  either  the  displacement 
of  the  center  of  mass  (except  for  a  rigid  inclusion)  or  the  shesr  stress 
at  the  center. 

Furthermore,  we  will  make  some  comments  on  the  difference  between 
the  small-time  asymptote  and  the  large-time  behavior  for  a  step  input,. 

When  ttie  incident  wave  encounters  the  Inclusion  at  t  “  -1,  the  inter¬ 
nally  transmitted  coapressional  wave  at  normal  incidence  to  the  surface 
will  reach  the  center  first  at  t  *  :  -  1,  reflect  at  the  leeward  surface, 
and  ratum  to  ths  center  at  t  -  3c  -  1.  The  expression  Eqs.  (25a)  or 
(25b)  essantially  only  treat*  t'\e  res*.  ?nsc  from  this  normal-incidence 


f ' 


Fig.  il  — Corjecfute  on  the  qualifo'ive  (esponje  of  cn  embedded  inclusion 
o  e  to  on  irvcidenr  compression*:!  step  wave 
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wave  and  neglects  the  shear-mlsmat ch-induced  waves  ( u  4  3.)  when  the  in¬ 
cident  wave  is  oblique  to  the  spherical  surface.  Thus  the  reflections 
in  Figs.  8,  9,  and  10  probably  do  not  give  the  total  response  after  one 
or  two  reflections,  except  if  c  <<  1.  If  we  may  be  permitted  to  con¬ 
jecture,  the  total  response  of  the  pressure  «t  the  cente-  of  the  inclu¬ 
sion  for  a  step  input  may  resemble  Fig.  1  A  for  c  n,  0(1).  The  small¬ 
time  overshoot  T  and  large-time  overshoot  may  be  useful  as  guides  to 
designing  °uch  a  pressure  tra^  aducer  or  interpreting  the  results  of  an 
actual  test. 

Finally,  we  should  mention  that  actuil  materials  have  Internal 
friction  or  damping  properties  which  were  neglected  here.  To  include 
these  properties,  the  elastic  moduli  can  be  replaced  by  the  complex 
counterparts  as  a  function  of  frequency  in  F.q .  (14).  Material  attenu¬ 
ation  of  sound  waves  would  be  Important  in  a  transducer  made  f  plastic, 
e.g.,  polyethylene  (Figs.  4  and  9),  and  in  fact  might  be  desirable. 
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VI.  ON  Thr.  DESIGN  OF  A  STRESS-STRAIN  TRANSDUCER 


At  the  beginning  of  this  Memorandum  ve  outlined  the  important 
characterisric3  that  a  sTess-strain  transducer  should  have.  In  the 
subsequent  sections  we  analyzed  t'ae  trans.er.t  response  of  a  spherical 
pressure  transducer  and  discovered  the  parameters  that  affect  the  r _  - 
sponse  characteristics.  We  propose  here  to  describe  one  possible  do- 
si  for  such  a  trs-scucer.  which  may  serve  as  a  point  of  reference 
for  all  apparent  areas  of  field  testing,  with  the  pc-seio1'*  exception 
of  ground-shock  testing. 

OVERALL  DESIGN  CONCEPT 

The  proposed  stress  transducer  will  sense  the  pre».  ;re  in  a  ground- 
shock  environment  by  measuring  the  d* '  dsMon,  or  a\  '"age  strain  at 
the  center  of  a  spherical  elastic  incision  h.  means  of  strain  gages. 

The  sphere  is  to  be  bonded  to  a  cylindrir'  .  sample  of  the  ground  roate- 
i  and  grouted  in  a  bore  hole  (Fig.  12).  Tr.e  spherical  shape  and 
choice  of  average  strain  will  result  in  a  transducer  insensitive  to 
orientation  relative  to  the  wave  front  ...  rvature  of  the  wave  front. 

This  unorthodox  she.  a  may  create  some  ratr. ctlon  problems.  Sev¬ 
eral  a i  .arnatives  are  availab:  •.  The  sphere  may  be  split  into  hemi¬ 
spheres,  the  gages  instal’ed,  and  the  sphere  fastened  or  cemented  to¬ 
gether.  Thi3,  however,  might  destroy  tu.  spherical  symmetry  of  its 
response.  If  the  inclusion  material  has  a  low  melting  point  or  can  be 
solidified  from  the  liquid  state  at  low  or  room  tomp  ratures  (e.g., 
epoxy),  the  gages  might  then  be  directly  embedded  in  the  sphere  before 
solidification.  This  method,  although  attractive,  will  usually  mean  a 
large  mismatch  in  acoustic  impedance  pc  between  inclusion  and  matrix, 
and  hence  large  internal  reflections  and  -  long  transient  decay  time. 
Under  certain  r ' rcumstances  a  large  impedance  mismatch  might  be  pre¬ 
ferred,  especially  in  high  pressure  shocks  where  the  pressure  in  the 
inclusion  would  be  attenuated  be' the  fracture  lim’t  by  t.,e  choice 
of  a  soft  material  such  as  plastic  or  lead. 

A  third  alternative,  the  one  preferred  here,  ia  to  drill  jmai1 
radial  access  holes  through  the  center  and  to  cement  semiconc  etor 
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Upper  cylindrical 
m  itrix  half 


Coaxial  coble  '  >r 
tv>c  iead-shie'Jed 
cables  i 
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Hemispherical 
transducer  seat 


Lower  cylindrical 
matrix  half 


Assembled  matrix -inclusion  system; 
sphere  and  cylinder  halves  cemented  toqether 
(e.g.  epoxy-type  bonding) 


I 


t 


Fig.  12 — Installation  details  of  the  cylindrical  matrix  -  spherical  transducer  assembly 
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1 

i 

> 

i  strain  gages  along  the  axis  of  the  cylindrical  hole  several  hole  diam- 

f  eters  away  from  tne  center,  thereby  avoiding  stress  concentration  ef- 

[  fects  due  to  the  tri-hole  intersection.  If  the  access-hole  diameter 

;  is  small  compared  with  the  inclusion  rauius,  diffraction  effects  should 

be  small.  Frequency  response,  however,  will  be  limited  by  the  acoustic 
impedance  ratio  and  inclusion  radius  rather  than  the  gage-hole  diameter. 

MATERIAL,  STRENGTH,  AND  SIZE  REQUIREMENTS 

To  minimize  the  spurious  transients  in  the  transducer  which  are 
due  to  internal  reflections,  the  inclusion  material  should  be  chosen 
so  that  its  acoustic  impedance  is  as  close  as  possible  to  that  of  the 
ground  matrix.  For  example,  in  a  granite  matrix  (Table  1)  aluminum 
voul^  be  the  best  choice  for  inclusion  m  terial  to  minimize  overshoot 
and  reflections. 

In  regard  to  the  required  strength  of  the  transducer,  it  should 
be  remembered  that  the  pressure  at  low  frequencies  in  the  inclusion  may 
be  higher  or  lower  than  that  in  the  free  field.  For  a  granite  matrix 
P/?q(w  -*  0)  *  1.4  for  a  steel  inclusion,  P/Pq  »  1.3  foi  titanium,  and 
P/Pq  ~  0.81  for  magnesium.  Hi  us ,  while  a  magnesium  inclusion  can  sup¬ 
port  a  matrix  pressure  24  percent  above  its  fracture  strength,  a  steel 
inclusion  can  support  a  matrix  measure  only  72  percent  of  its  fracture 
strength , 

The  size  of  the  device  is  .limited  by  the  seismic  bore  hole,  the 
size  of  the  strain-gage  access  holes,  and  the  desired  frequency  response. 

There  are  two  "figures  of  merit"  for  the  transient  response.  One 
is  the  rise  time  or  time  to  transient  decay  due  to  a  step  input,  the 
other  is  the  frequency  response  as  determined  from  the  pressure  ratio 
versus  wave  number  or  frequency  curves  (Figs.  2  through  5).  The  rise 
time  or  transient  damping  time  is  Independent  of  the  si^u  and  depends 
only  on  the  acoustic  impedance  ratio  of  matrix  and  inclusion. 

The  frequency  response  is  usually  defined  by  the  frequency  at  which 
the  gain  A(u)  5  P/PQ  is  up  or  down  3  dB  (dB  =  20  log^A)  or  where  A  ■ 
1.41  or  0.71.  This  frequency  will  depend  on  the  radius  of  the  inclu¬ 
sion.  For  a  lead  inclusion  (4-in.  diameter)  in  a  granite  matrix,  the 


gain  or  pressure  ratio  ?/Pq  will  be  up  3  dB  when  the  frequency  is 
around  10,000  cycles/sec.  For  a  titanium  inclusion  the  frequency  re¬ 
sponse  is  around  0-47  kc/sec  for  a  4-in.  diameter  inclusion.  This 
3-dB  figure  may  be  a  little  misleading,  since  the  whole  response  for 
titanium  never  varies  more  than  about  4  dB  (Fig.  5),  and  for  aluminum 
(upper  curve  Fig,  2),  never  more  than  3  dB.  It  is  far  easier  to  extend 
the  frequency  response  by  matching  the  acoustic  impedance  than  by  making 
the  inclusion  smaller. 

STRAIN  GAGES  AND  ELECTROMAGNETIC  INTERFERENCE 

The  strain  sensors  for  this  transducer  are  semiconductor  strain 

(12) 

gages.  These  p is ze resistive  cylindrical  elements  would  be  cemented 
in  pairs  along  each  of  the  three  mutually  perpendicular  axes  of  the  ac¬ 
cess  holes  (Fig.  13).  The  gages  would  be  mounted  on  each  side  of  the 
center  of  the  sphere  at  a  distance  of  four  hole  diameters  away  from  the 
center  to  avoid  stress  concentration  effects  near  the  center.  Each  pair 
would  be  connected  in  series,  and  the  three  pairs  then  series-connected 
into  a  total  of  six  gages  forming  one  arm  of  c  Wheatstone  bridge.  The 
output  of  the  bridge  would  be  proportional  to  twice  the  dilatation; 
hence  the  pressure  at  the  center  of  the  sphere. 

If  transient  electromagnetic  radiation  is  a  problem,  the  gages  may 
be  isolated  from  the  sphere  which,  if  a  conductor,  will  check  the  in¬ 
trusion  of  electromagnetic  radiation  into  the  interior  tn.cugh  the  skin 
effect.  All  but  quasistatic  fields  would  be  shielded  from  the  strain 
gages. 

FABRICATION  AND  INSTALLATION 

The  smalj  size  of  the  gage-access  holes  will,  of  course,  present 
some  fabrication  problems.  However,  the  use  of  specially  designed 
tweezers  and  probes  and  a  low-powered  microscope  should  make  the  prob¬ 
lem  tractable.  Small-diameter  enamel-coated  lead  wire  should  be  threaded 
through  the  holes  before  installing  the  gages.  The  gages  can  be  connec¬ 
ted  ,n  pairs  and  installed  as  such,  so  that  all  the  leads  exit  from 
only  one  access  hole.  With  the  leads  soldered  as  one  arm  of  a  bridge. 
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only  two  leads  need  be  carried  back  to  a  ground  station  into  a  bridge 
through  a  coaxial  or  shielded  cable.  If  electromagnetic  radiation  is 
not  a  problem,  one  of  the  leads  may  be  grounded  to  the  sphere  itself. 

An  exploded  diagram  of  the  installed  transducer  is  shown  in  Fig. 

12.  Rather  than  simply  grouting  the  sphere  in  place  with  loose  ground 
material,  we  suggest  a  cylindrical  sample  of  the  ground  be  obtained  if 
possible  and  split  into  two  smaller  cylinders.  On  one  of  the  ends  of 
each  cylinder  a  hemispheric  hole  would  be  machined  to  fit  the  sphere. 
Through  one  of  the  cylinders,  a  small  axial  hole  would  be  needed  for 
the  transducer  cable.  The  two  halves  would  then  be  cemented  (with  some 
epoxy,  for  example)  to  the  sphere,  the  trarssducersphere  in  the  cylin¬ 
drical  matrix  would  then  be  lowered  into  the  seismic  bore  hole  and  grouted. 

There  are  of  course  many  other  factors  to  consider  in  properly  de¬ 
signing  an  instrumentation  system.  Some  of  these,  such  as  the  auxiliary 
electronic  en*'ipraent,  are  beyond  the  scope  of  this  study  and  the  expe¬ 
rience  of  its  authors.  However,  no  proper  evaluation  of  the  problems 
can  begin  without  a  specific  device  at  hand  with  concrete  sizes  and 
specifications;  we  have  attempted  to  provide  these  number*,  and  proper¬ 
ties  for  those  with  intimate  knowledge  and  expert  nee  in  field  testing. 

The  overall  advantage  of  this  design  is  not  lower  cost  or  ease  of 
installation  (though  that  may  be  true),  but  that  its  shape,  internal 
construction,  and  installation  allow  the  response  and  field  results  to 
be  directly  compared  with  a  theoretical  analysis.  The  output  of  the 
transducer  cm  be  systematically  decoded  to  determine  the  pressure  in 
the  incident  wave  field. 
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